Magnetic hysteresis effects in superconducting coplanar microwave resonators 
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We performed transmission spectroscopy experiments on coplanar half wavelength niobium res- 
onators at a temperature T = 4.2 K. We observe not only a strong dependence of the quality factor 
Q and the resonance frequency / rcs on an externally applied magnetic field but also on the mag- 
netic history of our resonators, i. e. on the spatial distribution of trapped Abrikosov vortices in the 
device. We find these results to be valid for a broad range of frequencies and angles between the 
resonator plane and the magnetic field direction as well as for resonators with and without antidots 
near the edges of the center conductor and the ground planes. In a detailed analysis we show, 
that characteristic features of the experimental data can only be reproduced in calculations, if a 
highly inhomogeneous rf-current density and a flux density gradient with maxima at the edges of 
the superconductor is assumed. We furthermore demonstrate, that the hysteretic behaviour of the 
resonator properties can be used to considerably reduce the vortex induced losses and to fine-tune 
the resonance frequency by the proper way of cycling to a desired magnetic field. 



PACS numbers: 74.25.Qt, 74.25.Wx, 84.40.Dc, 03.67.Lx 



I. INTRODUCTION 

When Charles P. Bean introduced his famous model for 
the magnetization of hard superconductors in 1962 [1, 2], 
he probably has not foreseen, that once related magnetic 
history effects may be of importance for circuit quan- 
tum electrodynamics [3—5], quantum information pro- 
cessing [6] or single particle detection [7]. However, in 
this manuscript we will demonstrate, that by taking ad- 
vantage of magnetic hysteresis it is possible to tune the 
properties of coplanar superconducting microwave cavi- 
ties as currently intensely investigated and used in the 
above mentioned research areas. For many applications, 
the quality factor Q of the resonator, which defines the 
photon lifetime in the cavity and the sharpness of the 
resonance, is demanded to be rather high. Thus there 
are many current efforts to identify and suppress the dif- 
ferent loss mechanisms [8—11]. Recently advanced hybrid 
systems have been proposed [12-16], which couple ultra- 
cold atoms, molecules or electrons to superconducting 
microwave cavities or investigate the combination of ar- 
tificial atoms based on superconductors with real atoms. 
The experimental realization of these proposals is ex- 
pected to add another loss channel to the superconduct- 
ing microwave cavities, as the magnetic fields required for 
trapping and manipulation of the atomic systems [17, 18] 
will lead to the presence of energy dissipating Abrikosov 
vortices [19]. 

Lately, there were different approaches to reduce the 
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vortex associated energy losses in special experimental 
situations. In some experiments the magnetic field can 
be applied parallel to the plane of the superconduct- 
ing film, which reduces the flux in typical coplanar res- 
onators by orders of magnitude. This approach was used 
in experiments with spin ensembles that were coupled 
to microwave photons in superconducting transmission 
line cavities [20, 21]. For experiments that require out- 
of-plane magnetic fields it was demonstrated, that losses 
due to vortices can effectively be reduced by trapping and 
pinning the flux lines either in a slot in the center of the 
resonator [22], or in antidots that are patterned at the 
resonator edges [23] or all over the superconducting chip 
[24] . Patterning with antidots is particularly suitable for 
zero field cooling experimental conditions, when the vor- 
tices enter the superconductor from the edges and form a 
flux density gradient, called the Bean critical state. Un- 
der zero field cooling conditions there is also a number 
of hysteresis effects, which on the mesoscopic scale are 
related to the spatial distribution of Abrikosov vortices 



In this paper we will demonstrate that the character- 
istic properties of superconducting resonators (with and 
without antidots), i.e. the quality factor Q, the loss fac- 
tor 1/Q and the resonance frequency / rcs , can be tuned 
by taking advantage of spatial vortex redistribution. In 
particular we will show, that for fixed values of a perpen- 
dicularly applied magnetic flux density B the resonator 
losses due to vortices, 1/Q V (B), can be reduced and the 
resonance frequency f rC s{B) can be tuned by magnetic 
history. In a more detailed analysis we find strong indi- 
cations, that essential features of the measured hysteresis 
effects can only be described by assuming a highly inho- 
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mogcncous rf-current distribution in combination with 
field penetration models for thin films as first mentioned 
by Norris[27] and later discussed by Brandt and Inden- 
bom (NBI model) [28]. 

The paper is organized as follows. After this introduc- 
tory we first describe the sample fabrication and char- 
acterization techniques in Sec. II. Then, in Sec. Ill we 
present and discuss our experimental data that show 
a hysteretic behaviour in the vortex associated energy 
losses and the resonance frequency in perpendicular mag- 
netic fields. In Sec. IV we develop a simple model to 
describe the dependence of the vortex associated losses 
on the rf-current and vortex distribution in the resonator 
and compare our measurements with numerical calcula- 
tions. In Sec. V we discuss the possibility of exploiting 
the hysteresis to improve and tune the properties of the 
resonator in a specific magnetic field. Hysteresis effects at 
higher modes of the resonator and for non-perpendicular 
orientations between resonator plane and magnetic field 
are presented and discussed in Sec. VI . Finally, Sec. VII 
concludes the paper. 



II. RESONATOR FABRICATION AND 
CHARACTERIZATION 

We fabricated half wavelength coplanar transmission 
line resonators with a designed resonance frequency of 
/res = 3.3 GHz, that are capacitively and symmetrically 
coupled to feed lines via 90 /im wide gaps at both ends. 
Due to the small coupling capacitance, the resonators 
are strongly undcrcoupled[29]. In a previous study we 
investigated the influence of different spatial antidot ar- 
rangements on the performance of the resonators in mag- 
netic fields and thus patterned several resonators with 
and without antidots[23]. In the present study we inves- 
tigate these samples focusing on the magnetic hysteresis 
and its impact on the resonator properties. Figure 1 (a) 
shows a sketch of the resonator layout together with op- 
tical images of the resonators without (b), with one row 
(c), and with three rows (d) of antidots localized at the 
edges of the center conductor and the ground planes. The 
resonators are correspondingly named Res 0, Res 1+ and 
Res 3+ with a + indicating that the feed lines except 
for the rf launch pads are also perforated with antidots, 
cf. Fig. 1 (a). The antidots have a radius of R = 1 /im 
and an antidot-antidot distance of D = 4 fxm. For fur- 
ther design considerations of the antidots regarding size 
and arrangement, see Ref. [23]. 

All structures were fabricated on a single 330 /im thick 
2 inch sapphire wafer (r-cut) by optical lithography and 
subsequent dc magnetron sputtering of a d = 300 nm 
thick Nb film. After the sputtering, the wafer was 
cut into chips of 12 x 4 mm 2 , each one containing a 
single resonator. Finally the surplus Nb was lifted-off 
with acetone. The transmission line has a characteristic 
impedance of Z sa 540, where the width of the cen- 
ter conductor is S — 50 /im and the gap to the ground 




FIG. 1: (Color online) (a) Layout of a 12 x 4 mm chip with a 
capacitively coupled 3.3 GHz transmission line resonator and 
optical images of resonators with (b) 0, (c) 1, and (d) 3 rows 
of antidots. The parts of the feedlines, which are perforated 
with antidots are marked with + (cf. Ref. [23]). 



plane is W = 30 /im. The Nb film has a critical tem- 
perature of T c ps 9K and a residual resistance ratio of 
i?(300K)/i?(10K) ps 3.6. 

Each chip was mounted in a small brass box and 
the transmission line was electrically connected to sub- 
Miniature A (SMA) striplinc connectors using Indium as 
contact material. All measurements were performed at 
T = 4.2 K in liquid Helium (Helium-gas in the angle- 
dependent measurements). A magnetic field perpendic- 
ular to the resonator plane could be applied with a pair 
of Hclmholtz coils. We estimate the flux density seen 
by the resonator to be one order of magnitude larger 
than the applied external field due to flux focusing ef- 
fects. As we are interested in the resonator properties at 
magnetic fields of some mT, no measures were taken to 
shield the samples from earth magnetic field. Therefore, 
in all zero field cooling experiments we performed, some 
residual field was still present, but much smaller than 
any field we applied (except in the angle- dependent mea- 
surements). To characterize the resonators we applied 
a microwave signal of -P a pp = — 20dBm and frequency 
/ to one of the feed lines and measured the frequency 
dependent magnitude of the transmitted power Ptrans(/) 
with a spectrum analyzer. No attenuators or amplifiers 
were used in the measurements. We estimate the effec- 
tive power at the resonator input to be about 5-10 dB 
lower than at the generator output, to which all values 
of -Papp refer. Similarly, the power at the input of the 
spectrum analyzer Ptrans is about 5-10 dB smaller than 
the power directly at the resonator output. 



III. HYSTERESIS EFFECTS 

Figure 2 (a) shows the frequency dependent trans- 
mitted power Ptrans(J) around the fundamental mode 
n = 1 of a resonator without antidots (Res 0) for dif- 
ferent values of applied magnetic field between B = 
and B = 4 mT. Measurements shown were performed di- 
rectly after zero field cooling of the sample. The step 
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size of the applied magnetic field between the different 
curves is AB = 0.48 mT starting from B = with an 
additional curve at B — 4 mT. For better visibility adja- 
cent curves are shifted by +2dBm, where Ptrans(4mT) 
is the unshifted reference. As has already been re- 
ported before, we find that with increasing magnetic field 
the resonance frequency shifts downwards and the reso- 
nance peak gets smaller and broader, indicating increas- 
ing losses[19, 22, 23]. When we reduce the magnetic field 
from B = 4mT back to B = 0, we find a strong hys- 
teresis in the resonance characteristics. Figure 2 (b) de- 
picts the corresponding spectra (for the same values of 
applied magnetic field as in (a)), showing that the orig- 
inal state of (a) is not restored at B = 0. Interestingly, 
dfres/dB changes sign at B w 3mT: After a first in- 
crease of the resonance frequency with decreasing P, / rcs 
decreases again for B — >• 0. 
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FIG. 2: (Color online) Measured transmitted power Ptrans vs. 
frequency / of resonator after zero field cooling, when (a) 
the applied magnetic field was swept from to 4mT and (b) 
back to zero. Adjacent curves are shifted by +2dBm for bet- 
ter visibility, with Ptrans (4 mT) being the unshifted reference. 
AP ~ 0.48 mT, for details see text. The applied microwave 
power was P apP = — 20dBm. 

We fitted the measured transmission spectra with a 
Lorentzian and extract the resonance frequency / res (P) 
and the full width at half maximum A/(P) to quantita- 
tively analyze the hysteresis in the resonator properties. 
Using / re s(P) and A/(P) we calculate the magnetic field 
dependent quality factor Q(B) = / res (P)/A/(P) and 
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FIG. 3: (Color online) (a) Vortex associated energy losses 
1/Q V (B) and (b) resonance frequency / rca (P) of the funda- 
mental mode n = I of resonator for a full cycle of P. 



the magnetic field dependent losses per cycle l/Q(B), 
respectively. In order to quantify the losses associated 
with the magnetic field, i.e. the presence of Abrikosov 
vortices, we define l/Q v {B) = l/Q(B) - 1/Q(0), cf. 
Refs. [19, 22, 23]. Figure 3 shows the (a) vortex as- 
sociated losses 1/Q V (B) and (b) the resonance frequency 
fres(B) for a full magnetic field cycle, that is B = OmT 
-4 B = 4mT ^ B = -3.96 mT ->• B = 0. Arrows 
indicate the sweep direction of the field. 

Both the vortex associated losses 1/Q V (B) and the res- 
onance frequency / r cs(P) show a pronounced hysteretic 
behaviour. Interestingly, immediately after the sweep 
direction is changed at B = 4mT, the losses decrease 
considerably and are significantly smaller than the losses 
for the same field values, when coming from the virgin 
state. In the following we indicate properties, that refer 
to the upsweep of B with f and to the downsweep with 
j,. To emphasize the very first upsweep from the virgin 
state we use f|\ Over a considerable range of applied 
magnetic field 1/Q\,(B) < l/Qj(P) with a minimum of 
the losses at B w 2 mT. At B w 1.4 mT the loss curves 
from upsweep and downsweep cross and for even smaller 
magnetic fields l/Ql(B) increases again. For negative 
magnetic field values we find a similar behaviour, where 
l/Ql(B) < l/Qi(B). When we repeatedly sweep the 
magnetic field between ±P max = 4mT, 1/Q V (B) fol- 
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lows a butterfly like curve and never returns to the virgin 
state. In view of Abrikosov vortices an obvious interpre- 
tation of this hysteresis is, that we do not get the sample 
vortex free again during the field cycle. 

The resonance frequency in Fig. 3 (b) shows a hys- 
teretic behaviour very similar to the losses but inverted 
regarding absolute values. Immediately after the sweep 
direction is inverted, / res increases and after reaching a 
local maximum it decreases again with further decreas- 
ing magnetic field. Note, the magnetic field values at 
the downsweep, where 1/Q V (B) reaches a minimum and 
/res(-B) reaches a maximum are not the same. This dis- 
crepancy was observed for all resonators. We believe, 
that 1/Q V (B) and f res (B) have their extrema at different 
field values, because the frequency shift due to a change 
in the kinetic inductance has not only contributions from 
the normal conducting vortex cores but also from the 
global and local screening currents, which have a slightly 
different distribution than the vortex cores. 

Remanent vortices in type-II superconductors have al- 
ready been observed and investigated in various studies 
before, cf. [25] and references therein, leading also to 
hysteresis effects in the microwave properties of super- 
conducting structures [30-32]. In the following we will 
show, that by taking a closer look at the measured hys- 
teresis curves one can gain a new insight into the un- 
derlying physics, as the particular shape of the curve is 
intimately related to the microwave current distribution 
as well as the vortex distribution in the resonator. 



IV. THE RESONATOR LOSS MODEL 

In this section we introduce a simplified model, that 
allows us to derive an approximate expression for the 
dependence of the vortex associated losses l/Q v on a 
spatially varying flux density B(x) as well as on the 
microwave current density j (x) in a superconducting 
coplanar resonator. For simplicity we only consider the 
center conductor of our coplanar line and ignore the 
ground planes. Figure 4 shows a sketch of a coplanar 
waveguide with a center conductor of width S and thick- 
ness d. In the discussion in this section the magnetic field 
is always applied in the y-direction and the microwave 
current density always points in the z-direction. 

We start with classical mechanics and treat a vortex in 
the superconducting strip as a massless point-like particle 
under the influence of a driving Lorentz force Jl = j &o 
with the sheet current density j = j'q sin(cji) and a fric- 
tion fp = rjv, which leads to the one-dimensional equa- 
tion of motion rjv = j rl &o- In this picture the amplitude 
of the vortex velocity v is directly proportional to the 
amplitude of the alternating driving force and we find 
the energy dissipated per cycle AE to obey the propor- 
tionality AE = J Q T f L vdt oc (f Q { ) 2 with T = 2tt/uj. If 
we have several vortices at different positions Xi with 
local current densities jo(xi) we must sum up over all 
AE(xi) from single vortices to get the overall dissipation 
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FIG. 4: (Color online) (a) Sketch of a coplanar waveguide 
with center conductor width S, ground-to-ground distance a 
and thickness d; (b) Microwave current density distribution 
fo{x) normalized to the homogeneous distribution jhom = 
I/S, where / is the total current on the center conductor of a 
coplanar waveguide for different values of \l/S according to 
Eqs. (2), (3). 

A^totaJ = Ei ^E{ Xi ) cx Jo f (^) 2 - 

The flux density B(x) (divided by the flux quantum 
<f>o = 2.07 x I0~ 15 Tm 2 ) can approximately be treated 
as a continuous representation of the vortex density in 
the superconductor and with this we can rewrite the to- 
tal dissipated energy per cycle and per unit length in 
y-direction as 

A£ total cx ( '' \B(x)\(f { (x)) 2 dx:=5e(B), (I) 

J-S/2 

where | reflects the independence of the dissipa- 

tion contribution from the vortex polarity. The herewith 
defined quantity 8e{B) contains all information on the 
spatial distributions of vortices and driving forces and is 
hence feasible for the qualitative study of the experimen- 
tally found hysteresis curves, although it is not the same 
as 1/Q V (B). However, within our approach, there is the 
direct proportionality 1/Q V (B) oc 6e(B) and the propor- 
tionality factors are eliminated below as we only consider 
the normalized quantity Se(B) / <5e max . 

Another possibility to find the same proportionalities 
between Se, j (x) and B(x) from a more electrotechni- 
cal point of view starts with the dissipation power density 
Pd{x) = Re[p„(x)]jQ (x) 2 , which must be integrated over 
the width S of the cross section to get the overall dissipa- 
tion per unit length D = J S ^ 2 pd(x)dx. The hereby used 

real part of the vortex resistivity p v (x) can be obtained 
from the formulae given by different authors in different 
models [36-38] . All of these models have the linear pro- 
portionality p v (x) oc B{x) in common, what again leads 



5 



us to the above expression for Se(B). This argument is 
closely related to the one, which can be found in [19] to 
obtain an expression for 1/Q V (B). Interestingly, the ex- 
pressions in these models are derived treating the vortices 
as a vortex crystallite and not as individual particles and 
by adding a pinning potential. However, these assump- 
tions do not change the proportionalities between l/Q v , 
f f (x) and B(x). 

According to Ref. [26] and references therein the mi- 
crowave current density distribution Jq{x) in the center 
conductor of a coplanar line can be approximated by 



Jo f (*) 



A" 



(2) 



where 



C(z) 



1- # 



1 - 



0<|-| 
\x\ < 



< A, 



eff ■ 



(3) 

Here, a = S + 2W denotes the distance between the 
ground planes; in our case a = 110 /im. K is the com- 
plete elliptic integral and I is the total current. A c ff = 
Al coth(d/AL) is the effective penetration depth and Al 
is the London penetration depth [39, 40]. For \ e f[/S <C 1, 
the current distribution (2) is very inhomogeneous and 
has pronounced maxima at the edges of the strip. How- 
ever, with increasing X e s/S the maxima continuously de- 
crease until Jq{x) becomes completely homogeneous, i.e. 
Jo(x) = I/S = const., for X c g > S/2. Using the param- 
eters of our samples we find A e ff ~ Al = 100 nm and an 
approximate ratio of effective penetration depth to width 
of the center conductor of \ c a/ S = 0.002. Figure 4 (b) 
shows the current distribution according to Eq. (2) for 
four different ratios X h /S = 0.002,0.006,0.02 and 0.06. 
Note, we chose S and d in correspondence to our sam- 
ple parameters and varied Al instead. Therefore, A c ft is 
different for each curve with X cS /S « 0.002,0.008,0.069 
and 0.602. 

To describe the magnetic flux density B(x) in the cen- 
ter conductor of the resonator we start with the classical 
Bean profile[l]. During an upsweep from the virgin state, 
the flux density decreases linearly from the edges of the 
strip and can be expressed as 



B^(x,B ext ) = 



'^\x\ - (B* -Bext), §>\x\>§b, 



0, 



\x\ < f b, 



(4) 

where b = (1 — ^jp 1 ) and B* represents the applied field, 
when the flux fronts from both edges of the strip meet at 
x = 0. After the virgin upsweep to -B max , the flux profile 
for the downsweep is given by 



£ £ m ax gcx t ) (5) 



The flux density profiles B^(x, B ext ) and B^-(x, B ext ) are 
schematically shown in Fig. 5 (a) and (b) for several val- 
ues of -Bext during (a) a field upsweep to -B max = B* 




-0.5 0.0 0.5 
B IB* 



0.5 0.0 0.5 

B IB* 



FIG. 5: (Color online) Classical Bean model flux density 
B/B* in a superconducting strip of width S during (a) the 
upsweep to B cxt = B* (B oxt /B* = 0.2,0.4,0.6,0.8 and 
1) and (b) during downsweep from _B cxt = B* to £> ex t = 
—B* (B ext /B* = 1,0.6,0.2,-0.2,-0.6 and -1); Calculated 
<5e(B)/<5e max during a magnetic field cycle B cx t/B* = — s> 
Bcxt/B* = 1 — > B cxt /B* — — 1 — > Bcxt/B* — assuming the 
classical Bean flux density and (c) a homogeneous as well as 
(d) a highly inhomogeneous microwave current density with 
\ L /S = 0.002 and d/S = 0.006. Arrows in (a, b) indicate 
the progression of the subsequent flux profiles, arrows in (c, 
d) indicate the sweep direction. 



and (b) during the downsweep from B cxt = B* to 
Bext = —B*. Note, for each applied magnetic field B ext , 
\B^(x)\ < | -B e xt ] and \B±(x)\ > \B^(x)\. The aforemen- 
tioned relation also reflects on the dissipation, such that 
for the same values of applied magnetic field the losses 
during the downsweep should be larger than in the up- 
sweep. 

The lower part of Fig. 5 shows the calculated quantity 
(5e/<5e max for the Bean model flux profile with a homoge- 
neous (X e g/S > 1, (c)) as well as a highly inhomogeneous 
(A e ff/S = 0.002, (d)) microwave current distribution. All 
calculations of de(B) were carried out numerically with 
a spatial resolution of Ax = 2 nm. We repeated our cal- 
culations with different spacings (lnm < Ax < 20 nm) 
and found relative deviations < 3%. 

In case of a homogeneous microwave current density, 
the dissipation only depends on the number of vortices in 
the sample, not on their spatial distribution, cf. Eq. (1). 
Starting in the virgin state, the total amount of flux 
increases quadratically with applied field. Therefore, 
Se^(B) also has a positive curvature. As expected for 
this field profile (see above), Se^ > 5e^ and the mini- 
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mum value of 5e^, which corresponds to the least amount 
of flux in the sample, is reached for a negative value of 
applied field. 

In case of a highly inhomogeneous current distribution 
[Fig. 5 (d)] only a small region near the edges of the center 
conductor is responsible for almost all of the dissipation, 
cf. Fig. 4(b). The flux density in this area is almost iden- 
tical to i? GX t, hence the hysteresis is much smaller than 
in the homogeneous case. Although the hysteresis loop 
is significantly smaller, certain characteristic features of 
Se(B) remain unchanged, such as the positive curvature 
of Se^(B), the relation Se^ > Se^ and the fact, that the 
position of the downsweep minimum is in the negative 
field range. Consequently, the classical Bean field profile 
does not lead to a hysteresis as observed in our experi- 
ments. 

An alternative model of the flux density in thin film 
geometries was first considered by Norris [27] and later 
discussed by Brandt and Indenbom [28]. It is basically 
the aforementioned Bean model adapted to the geometry 
of thin superconducting strips. In this Norris-Brandt- 
Indenbom (NBI) model the flux density is given by 

{i-l y/(x-S')(x+S') 5 I i „, 
taim Mta„h(iWB c )> 2 > \ X \ > b 
0, \x\ < S' 

(6) 

where S' — S/2 cosh(B cxt / B c ), B c = fiodj c /n is the char- 
acteristic field, cf. Rcf. [28] and j c is the critical cur- 
rent density of the superconductor. Using the estimated 
j c w 5 x 10 10 A/cm 2 of our Nb, we calculate aB c «6 mT. 
As can be seen from Eq. (6), the NBI model leads to an 
excess flux density at the edges of the strip compared to 
the classical Bean profile. The downsweep flux density 
B^(x, B cxt ) is again defined according to Eq. (5). 

The NBI flux density profile is shown in Fig. 6(a) for 
several values of B cxt during a field upsweep to B cxt = 
3B C and (b) during the downsweep from B ext = 3B C to 
B C xt = SB C . As can be seen in Fig. 6 (b), B^ shows a 
remarkable behavior. For an externally applied field that 
is still positive, the flux density close to the edge of the 
thin film is already zero. When B cxt is decreased further, 
B^ at the film edge becomes negative and the point of 
zero flux density moves deeper into the sample, separat- 
ing areas with anti- vortices from that with vortices. 

To achieve comparability with experimental results, 
where B max = 4mT, a field range of |-B C xt| 1^B c k, 6mT 
was chosen for numerical calculations. Figures 6 (c,d) 
show Se(B)/Se niax for a (c) homogeneous and a (d) 
highly inhomogeneous current density. In the calcula- 
tions we avoided the divergence of the NBI flux density 
at the strip edges by positioning them between two in- 
tegration points, i.e. by effectively introducing a cutoff 
for B at Ax/2 from the conductor edges. In case of 
Jo( x ) = const., the NBI and the classical Bean model 
lead to similar Se(B)/Se max dependences, cf. Fig. 5 (c), 
which — as already mentioned — disagree with our ex- 
perimental data. For the inhomogeneous rf current dis- 
tribution, however, 5e(B)/Se max reproduces almost all 
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FIG. 6: (Color online) NBI model flux density B/B c in a 
superconducting strip of width S during (a) the upsweep to 
Bcxt = W c (Boxt/B c = 0.6, 1.2, 1.8, 2.4 and 3) and (b) during 
downsweep from _B oxt = 3B C to B cxt = — 3B C (B cxt /B c = 
3,1.8,0.6,-0.6,-1.8 and -3); Calculated <5e(B)/<5e max dur- 
ing a magnetic field cycle B ex t/B c — — > B cxt /B c — 1 — > 
B ex ±/B c = — 1 — > B cx t/B c = assuming the NBI flux den- 
sity and (c) a homogeneous as well as (d) a highly inhomo- 
geneous microwave current density with Xl/S = 0.002 and 
d/S = 0.006. Arrows in (a, b) indicate the progression of the 
flux profile, arrows in (c, d) indicate the sweep direction. 



characteristic features of the measured curve, cf. Figs. 3 
(a) and 6 (d). In particular, the hysteresis loop has a 
butterfly like shape, where Se^(B) < Se^(B) for a consid- 
erable range of B oxt . Also, B" exhibits a predominantly 
negative curvature and the minimum of 8e^(B) can be 
found at B cxt > 0. 

Yet, there is also a difference between experiment and 
theory. An asymmetry between the first very slow in- 
crease of the losses at small fields and the abrupt decrease 
immediately after the inversion of the sweep direction is 
clearly visible in experiment, but not in theory. Using the 
models described before, the curvatures d 2 B^(x)/dB^ xt 
at -B ext = and d 2 B ir (x)/dB c 1 xt at i? cx t = B max are di- 
rectly linked, cf. Eq. (5). Therefore, there is always a 
symmetry between the first slow increase of the losses 
at the beginning of the upsweep and the slow decrease 
at the beginning of the downsweep. We believe that the 
asymmetry in experiment originates from the small but 
nonzero lower critical field B c \ of the superconductor, 
which is not included in the theory. The existence of B Q \ 
inherently leads to an asymmetry between the upsweep 
from the virgin state and any following sweep. 

In Fig. 7 the dependence of the hysteresis loop on 
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FIG. 7: (Color online) (a) 5e(6ffB e xt)/^e mB x vs external 
flux density B ext /B c calculated for the NBI model with 
d/S — 0.006 and different flux focusing factors bff. (b) 
<5e(6fi_B ex t)/(5e m ax vs. cycled flux density B ext /B c calculated 
for the NBI model with d/S — 0.006 and different ratios Al/S 1 . 
Adjacent curves are subsequently shifted by +0.6 for better 
visibility. 



the ratio Al/S" and the flux focusing factor bg = 
B c g/ B cxt with the effective flux density B c g in the center- 
conductor-to-groundplane gaps is shown. Figure 7 (a) 
shows 5e(bffB ext )/8e ma _ x for four different flux focusing 
factors bg = 1,2,5 and 10. All curves exhibit the same es- 
sential features of the hysteresis but with some minor dif- 
ferences. With increasing flux focusing, the field regions 
around B cxt = and B cxt = B max where curvatures of 
5e^ and Se^ change sign become effectively compressed. 
Furthermore, the position of the downsweep minimum 
shifts with bg. By comparing calculations with exper- 
imental data, we find best qualitative agreement for a 
flux focusing factor between bg pa 2 and bg pa 5. Al- 
though the estimated value of B c pa 6 mT is comparable 
to -B m ax = 4mT in experiment and a flux focusing fac- 
tor of bg pa 5 is compatible with the resonator geometry, 
we would like to emphasize, that the theoretical model 
presented in this paper is too simple to be used for a 
quantitative analysis of our experimental data. 

Besides the flux focusing factor the homogeneity of 
the rf current distribution described by X L /S deter- 
mines the shape of the hysteresis loop. Figure 7 (b) 
shows Se(B)/Se niax for four different ratios Xl/S — 



0.002,0.006,0.02 and 0.06, the corresponding current 
distributions are shown in Fig. 4 (b). Starting from 
Xl/S = 0.002, which corresponds to our experimental 
conditions, the hysteresis loop becomes smaller with in- 
creasing homogeneity. For Xl/S sa 0.02 hardly any hys- 
teresis can be seen. When Al is increased further, up- 
and downsweep curves do not even cross anymore. 

Due to the good agreement between theory and exper- 
imental data, some important implications of our model 
and Se(B) calculations presented in this paper should bc 
addressed. In our calculations we have assumed, that 
the microwave current density and the flux (vortex den- 
sity), can be treated independently from each other. In 
general, however, the presence of a Bean-like flux gradi- 
ent leads to a redistribution of a transport current and 
vice versa [28]. Also, the oscillation amplitude of the 
vortices Ax is assumed to bc much smaller than the 
length scale on which the current distribution varies, such 
that Jq{x ± Ax) sa f${x). We believe, that both as- 
sumptions are reasonable, as in experiment we observe, 
that the vortex associated losses 1/Q V are almost in- 
dependent of the power of applied microwaves, as long 
as -Papp < OdBm [23]. If the microwave self- field sig- 
nificantly disturbed and rearranged the static magnetic 
field configuration by introducing additional vortices, we 
would expect some nonlinearitics [42], which would lead 
to a power dependence of the losses. 

Further simplifications arc, that in experiment, the flux 
focusing factor depends on the applied flux density, as 
the external flux is only partially focused into the gaps 
of the coplanar waveguide. Some of the flux also pene- 
trates the superconductor in form of Abrikosov vortices. 
Moreover, due to the geometry of our samples with a 
meandered resonator line, bg is expected to vary along 
the resonator. Consequently, also the losses Se — Se(z) 
depend on position along the resonator. Finally we ig- 
nore the ground planes and that they might experience a 
different effective B at their edges and hence a different 
hysteresis as the center conductor. In the experiments 
all these effects are merged together, but we think that 
none of them fundamentally changes the observed and 
analyzed hysteretic behavior. 



V. DEMAGNETIZATION AND TUN ABILITY 

In order to explore the possibility to return to the 
virgin state after magnetic cycling, we performed a de- 
magnetization procedure, i.e. we swept the magnetic 
field up and down again with decreasing amplitude and 
monitored the corresponding resonances. The procedure 
slightly differs from commonly used demagnetization cy- 
cles, as we kept the step width of the magnetic field 
sweep constant and thus with decreasing amplitude of 
the sweep, the oscillation period between positive and 
negative maximum became smaller accordingly. Still, in 
the following we will refer to this procedure as demagne- 
tization. 
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FIG. 8: (Color online) (a) Measured vortex associated energy 
losses 1/Q V {B) of Res of the fundamental mode n = 1 for 
a full demagnetization cycle (line), for the virgin field sweep 
(full squares) and minimum values for positive magnetic fields 
(open squares); (b) 1/Q V (B) of the resonators 3+ (triangles), 
1+ (circles, shifted by +5 x 10 _J ) and (squares, shifted 
by +1 x 10 -4 ) for the virgin field sweep (full symbols) and 
the minimum values for each magnetic field values during the 
demagnetization sweep (open symbols). 



Figure 8 (a) shows the measured energy losses 1/Q V (B) 
during a full demagnetization procedure (blue line). As 
can be seen, with decreasing amplitude of the up- and 
downsweeps, the losses at zero magnetic field also de- 
crease and approach, yet do not fully return, to the vir- 
gin state value. At the end of the demagnetization a 
small number of vortices remains in the resonator. Dur- 
ing each second half cycle Abrikosov anti- vortices ( "anti" 
relative to the remanent ones from the first half cycle) are 
pushed into the sample almost as far as the remanent vor- 
tex front reaches. The remaining losses after demagne- 
tization probably reflect remanent vortices (and possibly 
anti-vortices as well) which were not annihilated during 
the repeated cycle. To avoid non-equilibrium effects and 
to ensure an adiabatic magnetic field sweep, we waited 
about one minute after each field step. 

The demagnetization curve in Fig. 8 (a) shows, that 
over almost the whole magnetic field range the enve- 
lope of the full curve, i.e. the minimum values of the 



FIG. 9: (Color online) (a) Measured resonance frequency 
fica(B) of Res for a full demagnetization cycle (line) with 
maximum values (full squares) and minimum values (open 
squares) for positive magnetic fields; (b) Tunability t rcs (B) = 
/ res , max(-B) - / res , min(-B) of the resonators 3+ (triangles), 
1+ (circles) and (squares). 



vortex associated losses (open squares) , lies considerably 
below the values of the virgin field sweep (full squares). 
This provides the opportunity to significantly reduce the 
losses of resonators or other microwave circuitry compo- 
nents operated in magnetic fields (e.g. when it comes to 
trapping of an ultracold atom cloud in the vicinity of a 
resonator) by proper choice of magnetic history. We are 
aware, that the achievable reduction is small compared to 
other approaches such as using antidots or slots, but nev- 
ertheless our procedure can be used to further improve 
the device performance by about 30%. 

The demagnetization cycle behavior of the resonance 
frequency / res , is depicted in Fig. 9 (a). The curve 
strongly modulates and covers a range of about 1 MHz 
for all magnetic field values. For a quantification, wc 
define the frequency tunability of the resonator as the 
difference between the maximum and the minimum reso- 
nance frequency t rcs = / ros , max - / reSl m in for each value 
of applied magnetic field. The resulting t ICS (B) of the 
three resonators with and without antidots is shown in 
Fig. 9 (b) for positive values of magnetic field. For all 
three resonators the tunability is about t res (B) = 1 MHz 



9 



and hence seems to be almost independent of a perfora- 
tion and magnetic field with a small tendency to increase 
with the number of antidots. 

A tunability like this might be useful, e.g. when 
it comes to a finetuning of the superconducting cavity 
to the (fixed) transition frequencies of ultracold atom 
clouds. As each frequency within the tunability range 
at each field point is accessible by at least two different 
histories and the maxima/minima of the frequency /loss 
hysteresis are not at the same field values, one has to 
check all possibilities to find the optimal combination of 
desired frequency and losses. As the exact parameters 
of the hysteresis slightly differ from device to device and 
probably from setup to setup, we can not give a common 
recipe for finding the best working point here, but each 
device has to be pre-characterized in the corresponding 
experimental situation. 

We emphasize, that we do not propose to tune the 
properties of planar superconducting microwave compo- 
nents by just applying a magnetic field and introducing 
Abrikosov vortices here. If one just would like to fab- 
ricate a frequency tunable device for instance, other ap- 
proaches with larger ranges and tunability velocities seem 
more promising [43-45]. But if the microwave compo- 
nents have to be operated in specific magnetic fields any- 
way, magnetic history effects provide a nice additional 
possibility to reduce the losses and tune the resonance 
frequencies by spatially re-arranging Abrikosov vortices. 



VI. HIGHER MODES AND OTHER ANGLES 

So far we have only considered the fundamental mode 
of our resonators and an operation in a perpendicular 
magnetic field. In the following we will present measure- 
ment results on higher modes and discuss the influence 
of the tilt angle between resonator plane and magnetic 
field on the hysteresis loop. 

First, in Fig. 10 we exemplarily show the measured 
energy losses 1/Q V (B) of resonator 3+ for (a) the fun- 
damental frequency n = 1 and (b) , (c) , (d) the first three 
harmonics n — 2, n — 3, n = 4 in perpendicular mag- 
netic fields. Clearly, the characteristic features of the 
hysteresis loop can also be seen for higher modes. 

Interestingly, with increasing frequency a fine structure 
within the hysteresis curves emerges, see e.g. Fig. 10 (d) 
at B w ±1.5 mT. This fine structure might be related to 
inhomogeneously distributed losses along the resonator. 
Depending on the resonator mode n, nodes and anti- 
nodes of the microwave current standing wave probe dif- 
ferent parts of the resonator, i.e. they contribute dif- 
ferently to the overall losses. As already pointed out, 
there is a variation of the flux focusing factor and thus 
the effective flux density along the resonator due to its 
geometry. The flux focusing is smaller in the meandered 
lines than in the straight part at the midpoint of the res- 
onator, as the adjacent gaps of the meander allow for 
partial sharing of flux. Consequently, for different flux 




B [mT] B [mT] 

FIG. 10: (Color online) Measured vortex associated energy 
loss 1/Q V (B) for the four lowest modes n = 1, n = 2, n = 3, 
n — 4 of Res 3+ for a full cycle of magnetic field. 

focusing factors the downsweep minima might occur at 
different applied magnetic fields and the superposition of 
the l/Q v curves of several antinodes (which is, what we 
measure) might lead to a fine structure on the hysteresis 
loop with a multiple minima. 

We also performed measurements, where the magnetic 
field was not applied perpendicular to the resonator plane 
but with an angle a ^ 90° (cf. Fig. 4 [a]). To do so, we 
used a different setup, where we rotated the sample in 
the field of a stationary superconducting high-field split 
coil. 

By rotating the sample and measuring the field de- 
pendent quality factor for different angles we determined 
the perpendicular and parallel orientations between mag- 
netic field and resonator, assuming that for a given mag- 
netic field losses are maximal for perpendicular and min- 
imal for parallel orientation. Note, a comparison of the 
measured losses for the perpendicular orientations in the 
two setups showed, that they are practically identical, 
although the zero field quality factor is smaller in the 
high-field setup probably due to a remanent magnetiza- 
tion of the cryostat components. 

Figure 11 shows the measured dependence of the qual- 
ity factor Q on the applied field B of the fundamen- 
tal mode of resonator 3+ for three different angles. In 
perpendicular fields, a = 90°, we find the well-known 
and already described hysteresis with significantly re- 
duced losses l/Q v (increased quality Q) on the down- 
sweep branch, see Fig. 11 (a). For a < 90°, the hystere- 
sis loop hardly changes, as can exemplarily be seen in 
Fig. 11 (b), where a ss 15°. Note, here the field range 
is about four times larger compared to the measurement 
with perpendicular orientation (a) and the general shape 
of the hysteresis and the downsweep improvement in Q 
are clearly visible. As sin(15°) ~ 0.26, measurement re- 
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FIG. 11: (Color online) Typical measured hysteresis curves 
of the quality factor Q of the fundamental mode n = 1 (res- 
onator 3+) during a full cycle of the magnetic field for three 
different angles (a) a — 90°, (b) a « 15° and (c) a ~ 0° 
between resonator plane and applied field. 



suits strongly suggest, that the resonator losses are pri- 
marily determined by the component of B perpendicular 
to the sample. 

The situation changes dramatically for a ~ 0°, where 
the magnetic field was swept between ±0.2 T, as can 
be seen in Fig. 11 (c). For magnetic flux densities 
B < 80mT the decrease in Q is relatively small and 
also reversible (not shown). When the magnetic field is 
increased further, the quality factor rapidly decreases. 
After sweep direction inversion Q remains low until B < 
50 mT, when the quality factor increases and almost com- 
pletely recovers to the original value of Q(0). For nega- 
tive fields we observe qualitatively the same behavior. It 
is important to note, that < Q^. 

Interestingly, a comparison of measurement results 
with the theoretical curves presented in Sec. IV shows 
best agreement with the prediction of the classical Bean 
model and not the NBI model, cf. Fig. 5. Although on 
the first view one indeed would expect the classical Bean 
model to be the adequate description for the flux density 
profile in this experimental situation (field parallel to a 
superconducting plane), one has to be careful with this 
interpretation for two reasons. First one- has to consider, 
that the film thickness d = 300 nm is about three times 
the penetration depth (A c ff ~ 100 nm), and therefore only 
1.5 times the diameter of one Abrikosov vortex. Hence 



the Bean model as kind of mean field theory might prob- 
ably need a modification to be suitable for calculating 
the resonator losses in this case. The step-like struc- 
tures in Q(B) (e.g. at B ss ±0.1 T) might reflect some 
kind of a discrete and non-smooth flux entry into and 
exit out of the superconductor. Moreover, it is likely, 
that the measured hysteresis is also partly the result of 
a small misalignment between the applied field and the 
resonator plane. At flux densities of B m 100 mT, an 
alignment error of only 0.5° already introduces a field of 
w 1 mT perpendicular to the film, which is large enough 
to significantly reduce Q, cf. Fig. 11 (a). It remains to be 
stated at the end of this discussion, that the hysteresis 
loop seems not to allow for reducing the microwave losses 
in the resonator, if the magnetic field is applied close to 
parallel to the superconducting film. 



VII. CONCLUSIONS 

We experimentally investigated the properties of 
3.3 GHz superconducting coplanar transmission line res- 
onators in magnetic fields. We particularly focused on 
resonators, that were zero field cooled to T = 4.2 K and 
afterwards exposed to magnetic fields cycled in the milli- 
Tesla range. We measured the resonance frequency / ros 
and quality factor Q of the resonators and found strong 
hysteresis effects, which on the mesoscopic scale are due 
to presence of Abrikosov vortices in the superconduct- 
ing film and their spatial redistribution during the field 
cycles. Using a simple model for the vortex associated 
resonator losses we show that different combinations of 
microwave current and flux density distributions lead to 
characteristically different hysteresis loops. We find best 
agreement between experiment and theory for a current 
distribution strongly peaked at the resonator edges and 
a modified Bean flux gradient for thin films, as described 
by Norris, Brandt and Indcnbom. We have also shown, 
that the hysteresis may be used to improve the resonator 
performance for fixed values of applied magnetic field by 
proper choice of magnetic history. Accordingly, the res- 
onance frequency can be tuned by about 1 MHz, i.e. for 
our resonators a few zero-field linewidths at liquid He- 
lium temperature. Both, the reduction of the losses and 
the tunability of the frequency are found to be qualita- 
tively and quantitatively independent of a pre-pattcrning 
of the resonators with antidots, which in a previous study 
has been shown to be feasible to reduce vortex associated 
losses. Furthermore, we show, that the hysteresis can also 
be found for higher modes n = 2,3,4 of the resonators 
and for angles between 90° and a few degrees between 
the field direction and the resonator plane. In the par- 
allel orientation, the hysteresis showed a very different 
behaviour with presumably no possibility for a resonator 
improvement with magnetic history. 

As for many experiments in circuit quantum electrody- 
namics the resonators are operated in the Millikclvin and 
single photon regime, the effects presented here have to 
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be investigated under these conditions in further studies, 
but there are definitely no obvious reasons, why the re- 
ported hysteresis effects should qualitatively change with 
decreasing temperature and power. Other parameters, 
which are expected to change the nature of the hystere- 
sis and have to be considered for possible applications are 
the geometric dimensions of the transmission line and the 
thickness of the films, which together with the magnetic 
penetration depth have a strong influence on current dis- 
tribution and flux density profile. 
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